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^ Abstract: We have operated a liquid-ai^gon large-electron-multiplier time-projection chamber 

^ (LAr LEM-TPC) with a large active area of 76 x 40 cm and a drift length of 60 cm. This setup 

_0 represents the largest chamber ever achieved with this novel detector concept. The chamber is 

^ equipped with an immersed built-in cryogenic Greinacher multi-stage high-voltage (HV) multi- 

^ plier, which, when subjected to an external AC HV of ~1 kVpp, statically chai^ges up to a voltage 

I ^ a factor of ^^30 higher inside the LAr vessel, creating a uniform drift field of ~0.5 kV/cm over the 

full drift length. This large LAr LEM-TPC was brought into successful operation in the double- 
j>. phase (liquid- vapor) operation mode and tested during a period of ~ 1 month, recording impressive 

O three-dimensional images of very high-quality from cosmic particles traversing or interacting in 

ly-j the sensitive volume. The double phase readout and HV systems achieved stable operation in cryo- 

genie conditions demonstrating their good characteristics, which particularly suit applications for 
next-generation giant-scale LAr-TPCs. 
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1. Introduction 

The design concept of the novel Liquid Argon Large Electron Multiplier Time Projection Cham- 
ber (LAi- LEM-TPC) and its promising performances have been demonstrated experimentally with 
a detector prototype having an active ai^ea of 10 x 10 cm^ []I], 0]. This new type of LAr-TPC, 
operating in the double-phase (liquid-vapor) mode in pure argon, is characterized by its charge am- 
plifying stage (LEM) and its projective anode charge readout system. Situated in the vapor phase 
on top of the active LAr volume, these latter provide an adjustable charge gain and two independent 
readout views, presently each with a pitch of 3 mm^ With Tq time information provided by the 

'Thanks to the amphfication stage which yields enough charge to be shared among several readout electrodes, smaller 
pitch sizes can be considered. 
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LAr scintillation, they provide a real-time three-dimensional (3D) track imaging with dE/dx infor- 
mation and the detector acts as a high resolution tracking-calorimeter. A very high signal-to-noise 
ratio can be reached in the LAr LEM-TPC thanks to the gaining stage. This significantly improves 
the event reconstruction quality with a lower energy deposition threshold and a better resolution 
per volumetric pixel (voxel) compared to a conventional single-phase LAr-TPC [0], such as e.g. 
ICARUS [ED, ArgoNeut [U or J-PARC T32 [0]. In addition the charge amplification compensates 
for potential loss of signal-to-noise due to the charge diffusion and attachment to electronegative 
impurities diluted in LAr, which both become more important as the drift length increases. It is 
therefore considered to be a promising technology for the next-generation giant-scale underground 
detectors for neutrino physics and proton decay searches (GLACIER) [0, EL ED and for direct Dark 
Matter searches with imaging [Q]. 

In this paper, we report on the first operation of a LAr LEM-TPC prototype having an active 
area of 76 x 40 cm^ and a drift length of 60 cm, which is the largest chamber ever realized with this 
novel detector concept. The detector was tested successfully in the double-phase operation mode 
in pure argon, recording events of cosmic tracks traversing or interacting in the active volume. 
The analysis of its tracking performance which requires full 3D reconstruction of a large number 
of recorded tracks is reported elsewhere mm . and the focus of this paper is put on the immersed 
high- voltage (HV) system for the drift field and cold operation. 

The HV system is based on the technique of a cryogenic Greinacher flTTll (also known as 
Cockcroft- Walton finp ) HV multiplier immersed in LAr directly inside the detector vessel. It is 
an innovative technique to generate a HV for creating a drift electric field in a LAr-TPC, with 
some advantages compared to other methods used up to now whereby the very HV is externally 
generated and brought into the detector via a HV feedthrough (for a comparison see Ref. [1T3I1 ). 
This technique particularly suits a large- to giant-scale detectors, where a very HV of ~ 100 kV to 
~1 MV would be needed for a drift field of ~1 kV/cm over a length of ~1 m to ~10 m. 

The main advantages of the Greinacher HV multiplier are: (1) all the HV parts are immersed 
in LAr which has a large dielectric strength (~1 MV/cm), (2) thus feedthroughs for very HV ai^e 
not needed, (3) the circuit itself can provide all the stages needed to create a uniform field, thus the 
system has no resistive load, (4) thus the power dissipation is virtually zero and (5) this allows a 
low frequency (e.g. 50 Hz) of the AC input signal which is fully outside of the bandwidth of the 
charge amplifiers used for this type of detectors. One main disadvantage is the difficultly of access 
in case of a failure. The potential impact of this access on the pure liquid argon volume can be 
large if the active components are directly immersed in it^. 

To study in detail the functionality of such a system under practical conditions incorporated in 
a LAr-TPC, we built a 30-stage Greinacher HV multiplier integrated to the large-area LAr LEM- 
TPC prototype. We earned out the test experiment at CERN from October to December 201 1 using 
cryogenic LAr apparatus of the ArDM Project (CERN RE18) linil . which is an experiment for a 
direct Dark Matter search exploiting a 1 -ton-scale double-phase argon TPC [Hj. 

In Section we describe the experimental setup, i.e. the LAr LEM-TPC, the cryogenic appa- 

^It is still possible to imagine that the Greinacher circuit is located in a separately evacuable volume in contact with 
the main pure liquid argon volume during normal operation, but which would be disconnected and separately emptied 
during maintenance, with minimal impact on the main liquid argon volume, thereby still avoiding the need for a HV 
feedthrough. 
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ratus and the HV system. Results of a series of measurements performed in air at room temperature 
and in LAr for characterizing the Greinacher HV multipher are presented in Section 0. In Section ^ 
we report on the test experiment, where the entire LAr LEM-TPC with the HV system was tested 
successfully in the double-phase operation mode in pure argon. Finally we discuss and conclude 
on the results in Section 0. 

2. Experimental setup 

2.1 Design of tlie large-area LAr LEM-TPC 
2.1.1 LAr LEM-TPC design overview 

Figure |I] illustrates the structure of the LAr LEM-TPC hanging inside the ArDM detector vessel. 
The field cage, which defines the active LAr target and drift volume, is constructed with four PCB 
side walls in a shape of rectangular box with dimensions 76 cm (L) x 40 cm (W) x 60 cm (H). On 
the inner surface of each PCB 31 horizontal field shaping strips, which hereafter are called field 
shapers (FS), are formed by etching at a constant pitch of 20 mm. The top (FSO) and the bottom 
(FS30) field shapers have a width of 9 mm while the ones in between have 18 mm, with a gap of 
2 mm between neighboring ones. Negative high voltages linearly decreasing from bottom to top 
ai^e distributed to the field shapers so that a vertical drift electric field is created uniformly over the 
full drift length of 60 cm. 



Figure 1. 3D CAD model of the LAr LEM-TPC prototype installed in the ArDM LAr vessel. The field 
cage has dimensions 76 cm (L) x 40 cm (W) x 60 cm (H). The LEM charge readout system is mounted 
on top of the field cage. Negative HVs linearly decreasing from bottom to top are generated by the built-in 
Greinacher HV multiplier and distributed to 3 1 field shapers formed by etching on the inner surface of the 
side wall PCBs. A vertical drift electric field thus is created uniformly over the full drift length of 60 cm. 
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On top of the field cage two horizontal extraction grids are mounted with a gap of 10 mm in 
between. The grid is a 0.15-mm-thick stainless-steel mesh, where square holes with a size 2.85 x 
2.85 mm^ are etched at a pitch of 3 mm all over the active area. The lower grid is positioned at 
the top face of the field cage. For the double-phase operation mode the LAr surface is adjusted at 
the middle of the two grids. The liquid level at each of the four comers can be monitored with a 
precision of ~0.5 mm with the aid of four capacitive level meters. Ionization electrons produced 
by ionizing particles are drifted upwards to the liquid surface. These electrons are extracted across 
the liquid-vapor interface into the gas argon (GAr) phase with the aid of a strong extraction field 
of typically 3-4 kV/cm between the two grids. They are then collected by the charge readout 
system incorporating the Large Electron Multiplier (LEM), which is described in Section [2.1.2| . 
The bottom face of the field cage is covered by the cathode grid that is a stainless-steel mesh of the 
same type as used for the extraction grids. The top (FSO) and the bottom (FS30) field shapers ai^e 
electrically coupled to the lower extraction grid and respectively, to the cathode. 

The HV for creating the drift field is generated using a built-in 30-stage Greinacher HV mul- 
tiplier, which is described in detail in Section P3| . The circuit is integrated in the design of one 
of the side-wall PCBs and the components are mounted directly on the outer surface of the PCB, 
as can be seen in Figure 0. As already mentioned the generator itself avoids the use of a voltage 
divider, since each multiplying stage provides a characteristic DC voltage. The various multiplying 
stages generate a monotonously increasing potential to supply the electrodes surrounding the drift 
volume. In our setup, the DC output of each of the 30 Greinacher stages is connected to each field 
shaper via a wire through the PCB. 




Figure 2. Picture of the LAr LEM-TPC with readout sandwich and drift field cage, which are about to be 
inserted into the ArDM vessel. The stack of blue capacitors on the front side wall shows the Greinacher HV 
multiplier. 

Under the field cage two photomultiplier tubes (PMTs) are installed to detect scintillation 
light produced by charged particles crossing the LAr target (primary scintillation, SI), as well as 
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the secondary scintillation light (S2) produced via electro-luminescence by the electrons extracted 
to the GAr phase I1T5I1 . We employed Hamamatsu R11065 cryogenic 3" PMTs having a high 
quantum efficiency of ~30%, specifically designed for the use in LAr [ITHl . The quartz window 
of the PMT was coated with a wavelength shifter, l,l,4,4-Tetraphenyl-l,3-butadiene (TPB) [HTIl . 
in a Paraloid™ B-72 polymer matrix to detect deep ultra violet photons of the argon scintillation 
(peaked around 128 nm) . The PMT signal is used primarily for triggering and for determination of 
the event Tq absolute time. To protect PMTs from possible discharges from the parts at HV, another 
mesh electrode, which is grounded to the detector vessel, is inserted between the cathode grid and 
the PMTs. The distance between the cathode and the protection grid was 9 cm and that between 
the protection grid and the PMTs was 1 cm. 

2.1.2 LEM charge readout system 

The charge readout system of the LAr LEM-TPC is a direct extrapolation of the successful 10 x 
10 cm^ prototype \^ to an active area as large as 76 x 40 cm^. It is characterized by (1) a single 
1-mm-thick LEM amplification stage and (2) a 2D readout anode realized on a single PCB. The 
system has a multi-stage structure consisting of two extraction grids, a LEM, a 2D readout anode 
and two signal collection planes as illustrated in Figure 0. The distances between the stages and 
typical configurations for the potentials and the inter-stage electric fields are summarized in Table |T]. 
In typical configurations a positive HV of ~ 1 kV is applied to the anode while the extraction grid 
in liquid is operated at a negative HV of ~ —7 kV, resulting in a virtual ground within the LEM 
plane. 



76 cm 




Figure 3. Structure of the charge readout system incorporating a LEM (Large Electron Multiplier) and a 2D 
anode. Right-front corner is cut in the illustration to present different stages in the multi-stage structure. 

The LEM is a macroscopic hole electron multiplier built with standai^d PCB techniques^. It is 

■'it was manufactured at ELTOS Circuit! Stampati Professionali, Italy (http : / /www . eltos . com/). 
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Table 1. LEM charge readout system configurations. 




Distance to the 


Typical operating 


Field to the stage 




stage above (cm) 


potential (kV) 


above (kV/cm) 


Signal collection plane 2 


— 





— 


Signal collection plane 1 


2 


+ 1 




2D anode 


1 


+1 




LEM (top electrode) 


0.2 


+0.5 


2.5 


LEM (bottom electrode) 


0.1 


-3 


35 


Extraction grid in GAr 


1 


-4 


1 


Extraction grid in LAr 


1 


-7 


3 


Cathode 


60 


-31 


0.4 



a 1-mm-thick FR4 plate, double-side cladded with a passivated copper layer. Of the order of half a 
million holes 500 jxra in diameter are CNC (Computer Numerical Control) drilled through the plate 
at a pitch of 800 jxra between the centers of adjacent holes. After the PCB has been manufactured 
and drilled, the copper layers are further chemically etched to create 40 /im wide rims around the 
holes. A HV of typically 3.5 kV is applied across the two faces creating a strong electric field of 
35 kV/cm in the LEM holes leading to multiplication of electrons by avalanches. 

The 2D anode plane has two orthogonal sets of readout strips on the bottom face providing 
two independent readout views. On top of the strips for one view, those for the other view are laid 
with a thin (50 /im) polyimide (Kapton) insulating bed underneath. For each view, copper strips 
are formed by etching at a pitch of 600 /^m covering nearly the full area of 76 x 40 cm^. The 
readout pitch of 3 mm is obtained by bridging five strips at one end. The strips on top are narrower 
(120 /xm) than those lying underneath (500 /xm) optimizing for an equal charge collection by the 
two views operating at the same potential. A more detailed description of the 2D anode can be 
found in dH]. The strips for both of the views are tilted by 45° with respect to the sides of the 
rectangular plate, to have symmetric conditions for both views. 

We use two signal collection planes, between which HV decoupling capacitors (270 pF) ai^e 
connected. Each readout channel is routed to a collective 32-channel connector to a fiat cable. A 
surge arrester protecting the readout electronics from discharges is mounted for each channel. The 
cables bring the signals through a custom-made feedthi^ough on the upper flange, connected to the 
readout electronics [ITOl placed outside of the vessel. 

The LEM charge readout system is assembled as a multi-layered "sandwich" unit with pre- 
cisely defined inter-stage distances and inter-alignment. The entire unit can then be mounted to or 
dismounted from a field cage, which may have different configurations, e.g. drift lengths. 

2.2 ArDM cryogenic system 

The double-phase operation of the LAr LEM-TPC requires a sophisticated cryogenic LAr apparatus 
which fulfills demanding conditions: (1) the LAr target must be maintained at a very high purity 
which corresponds to an oxygen-equivalent impurity concentration less than 1 ppb, (2) the pure 
argon volume must be pumped down to a high vacuum (~ 10^^ mbar) before filling with LAr 
to reduce contamination of the argon due to outgassing, (3) filling through a purification filter is 
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obligatory, (4) to maintain and improve the purity during the test period (i.e. months), LAr must 
continuously be cleaned with the aid of a recirculation/purification system, (5) thermodynamical 
conditions, i.e. temperatures and pressures must continuously be controlled and (6) the system must 
be protected from possible hazards due to a large amount of cryogenic LAr in case of unexpected 
incidents. Therefore, we carried out the test experiment fully exploiting the cryogenic apparatus of 
the ArDM Experiment, which we had previously built and operated at CERN 

Figure presents a rendered CAD picture of the apparatus. The ArDM cryogenic system con- 
sists mainly of the detector vessel and the LAr recirculation/purification system. The two sections 
are connected with horizontal pipings at the top and the bottom. 




irculation/ 
ition system 



Figure 4. Overview of the cryogenic system of ArDM. 



The detector vessel has an inner diameter of 1 m and a height of 2 m containing approximately 
2 tons of LAr when it is full. In this test experiment the vessel was filled up to approximately 
two thirds of the full volume. The pure LAr volume, which is constructed as a ultra-high vacuum 
tight closed volume, is surrounded by a separate LAr volume serving as a cooling bath. The 
temperature of LAr and consequently the pressure in the cooling bath are maintained with the aid 
of the cryocoolers installed in the recondenser sitting on top of the recirculation/purification section. 
The cooling bath then is surrounded entirely by an insulation vacuum. The pure LAr is recirculated 
through a custom-made purification cartridge filled with activated copper powder with the aid of 
a cryogenic bellows LAr pump installed below the cartridge, and is purified continuously, i.e. the 
residual oxygen contamination is removed through a reaction O2 + 2Cu — )• 2CuO. For initial filling 
of the vessel with LAr another custom-made cartridge of a similar type with a molecular sieve and 
an active copper powder layer was used. In addition a gas recirculation system incorporating a 
commercial SAES getter^ is implemented in the system. 

The thermodynamical conditions, i.e. pressure and temperature, at different parts of the system 
ai^e monitored and logged continuously by a PLC (Programable Logic Controller) system. Further- 

^SAES Getters S.p.A.. http : / / www . saesgetters . com/ 
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more the PLC system allows manipulation of various active components such as vacuum pumps, 
electro valves, the bellows pump and heaters to regulate the cooling power of the cryocoolers. Al- 
though the cryogenic system has a number of rupture disks as an ultimate safety measure the PLC 
system provides early alarms and a range of interlocks to stabilize the conditions before reaching 
the rupture pressure. In addition, the system can be accessed via internet allowing a remote control. 

During the test experiment, the LAr temperature was maintained at 86.3 K. We reached a 
pressure stability better than ±1 mbar. With the electron drift velocity of 1.4 mm//is at 86.3 K and 
the drift electric field of 0.4 kV/cm it takes 430 /is for electrons to travel over the full drift length 
of 600 mm of the tested LAr LEM-TPC. We ran the bellows pump for two weeks at full speed 
and the oxygen-equivalent impurity concentration reduced from ~1.5 to ~0.5 ppb. The lifetime 
of drifting electrons (t) can approximately be calculated as a function of the oxygen-equivalent 
impurity concentration (p) [ITEIl : T [jUs] ^ 300/p [ppb]. The impurity concentration of 0.5 ppb 
leads to T 600 /is. 

2.3 HV system 

2.3.1 Cryogenic Greinacher HV multiplier 

The 30-stage cryogenic Greinacher HV multiplier was designed fully exploiting experiences gained 
in R&D of the ArDM HV system [1T3I1 . which is based on the same concept and has 210 stages. For 
the main circuit components, namely capacitors and diodes, ones of the same types as in ArDM 
were chosen because of their known reliability in LAr. Figure § illustrates the diagram of the 
circuit. Each capacitor symbol corresponds to one single polypropylene capacitor having 82 nF 
and each diode symbol to a chain of three avalanche diodes connected in series. The direction of 
the diodes determines the polarity of the DC output voltage which is negative in this case. A single 
stage is designed to hold as high as 2 kVDC, i.e. the circuit is able to output a total of —60 kVDC. 




Figure 5. Diagram of a 30-stage negative Greinacher circuit. The stage No. 3 is colored in blue to visualize 
the structure of a single stage. The DC voltage source, the capacitors and the diodes shown in red are used 
to shift the potential of the whole circuit by a certain value with respect to the ground. 



The output voltage is proportional to the amplitude of sinusoidal input voltage. Fully charged 
the voltage across each stage reaches ideally the peak-to-peak value of the input voltage (Vpp) and 
hence the potential at the ith stage equals to Vj = i • V^^. The DC output of the /th stage which is 
denoted as Vj in the diagram is connected to the field shaper FS/ (/ = 0...30). The potential of the 
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whole circuit can be shifted by adding a DC voltage source and several appropriate components, 
as shown in red in Figure ^ This is used to adjust the potential of FSO and of the extraction grid in 
LAr with respect to that of the charge readout system. 

The Greinacher circuit is integrated directly on the outer surface of one of the side-wall PCBs 
as can be seen in Figure 0. The shifting part is built on a separate PCB and is mounted near the 
Greinacher circuit. For the last 12 stages a HV resistor having 200 MQ. is inserted between the 
Greinacher stage and the field shaper to limit the current flowing out of the large capacitances in 
case of an unwanted discharge from the cathode or the field shapers at very HVs. For the double- 
phase operation mode the whole circuit is immersed in LAr. 

Table summarizes the parameters of the Greinacher circuit as described above. The specifi- 
cations of its components can be found in Table 0. 



Table 2. Specifications of the Greinacher circuit. 



Number of Greinacher stages , N 


30 




Capacitance per capacitor, C 


82 


nF 


Number of capacitors 


60 




Combined capacitance 


5.5 


nF 


Maximum rated voltage per stage 


2000 


V 


Maximum output voltage (Knax) 


60 


kV 


Stored charge at Vmax 


10 


mC 


Stored energy at Vmax 


10 


J 



2.3.2 AC power supply unit 

A dedicated power supply unit was developed to supply appropriate AC voltages for charging 
Greinacher circuits for LAr-TPCs. Its block diagram is shown in Figure ^ and the relevant specifi- 
cations are summarized in Table ^ 

The AC power supply unit is controlled and monitored through a USB interface by a com- 
mercial PC. The amplitude of the sinusoidal output voltage can be ramped up from to 2.7 kVpp 
(peak-to-peak) with a maximum current of 10 mA (r.m.s.). A fixed frequency of 50 Hz for the 
output signal is generated by a quartz oscillator. The ramp-up rate of the output voltage can be set 
in the range from 10 to 1000 V/s. The output voltage and current ai^e measured by a peak-to-peak 
detector and an absolute peak detector, respectively, and are read out via USB. On the front panel 
there is an analog output for each of the output voltage and current to allow monitoring of their 
waveforms on an external oscilloscope. 

The internal control unit can be enabled/disabled either by the USB interface or by a manual 
push-button switch on the front panel. The peak output current is set by the interface through the 
control unit. If a sudden discharge of the Greinacher capacitors generates an over-current in the 
secondary coil of the transformer, the control unit turns off the system. This current limit is set 
by a potentiometer on the front panel. The power amplifier is also protected for over-voltages. 
The system has never been damaged even though a few discharges occurred while testing several 
different Greinacher circuits. 
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Table 3. Components of the Greinacher circuit. 



Capacitor 

Double metallized film pulse capacitor, polypropylene dielectric 



Manufacturer 


Evox Rifa" 




Type 


PHE450 




Winding construction 


Two Section 




Lead spacing 


22.5 


mm 


Dimensions 


11.0 X 21.5 X 26.0 


mm^ 


Rated voltage DC 


2000 


VDC 


Capacitance 


82 


nF 


Diode 






High- voltage soft-recovery rectifier (avalanche diode) 




Manufacturer 


Philips Semiconductors^ 




Type 


BY505 




Package 


Rugged glass package 




Maximum ratings 






Peak reverse voltage (-65 ... +120 °C) 


2200 


V 


Non-repetitive peak forward current 


5 


A 


Repetitive peak forward current 


800 


mA 


Average forward current'' 


85 


mA 


Electrical characteristics 






Diode capacitance (0 V) 


2 


pF 



"Evox Rifa GmbH, http : / /www . evoxrif a . f i/ 
^NXP Semiconductors, http : / / www . nxp . com/ 
•^Averaged over any 20 ms period at 25 °C. 



A software to control the AC power supply unit was built based on LabView^. It has a graphical 
interface as shown in Figure [7], where one can type in the set value for the peak-to-peak output 
voltage and can start/stop charging by pressing a button. The ramp-up rate of the output voltage 
can also be set using the interface. The monitored value of the peak-to-peak voltage and a waveform 
of the output current for one cycle of the sinusoid are displayed on the window. 

2.3.3 Discharging system 

A similar system as reported in [1T3I1 was built to actively and safely discharge the Greinacher 
circuit. Its circuit diagram and the mechanics are presented in Figure 0. A mechanical switch 
realized with a rotating discharging rod steered by a rotary motion feedthrough makes contact 
between the cathode and a resistor chain of total of 1 GQ. which is connected to the top flange of 
the dewar at the other end. The resistance is chosen to discharge the circuit completely in the order 
of 10 s and to keep the peak discharging current small enough compared to the maximum rating for 

^LabView, National Instruments Corporation, http : / /www . ni . com/ labview/ 
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Figure 6. Block diagram of the AC power supply unit for charging the Greinacher circuit. 



Table 4. Specifications of the AC power supply unit. 



Output 

Frequency 
Voltage 

Maximum peak current (r.m.s.) 
Transformer 

Manufacturer 
Type 

Power rating 
Primary circuit (input) 

Maximum voltage (r.m.s.) 

Maximum current (r.m.s.) 
Secondary circuit (output) 

Maximum voltage (r.m.s.) 

Maximum current (r.m.s.) 
Monitoring output 
Voltage monitor 
Current monitor 



50 
... 2700 
10 

Josef Betschart AG" 
TET 84a/10852 
20 

30 
780 

2 
10 

1 

0.5 



Hz 
V 

Vpp 

mA 



VA 

V 

mA 

kV 

mA 

V/kV 
V/mA 



"Josef Betschart AG. http : / /www . betschart -trafo.ch/ 



the forward current through the diode. The cathode voltage and the discharging current decreases 
with an exponential decay at a decay constant T = 6.6 s which is determined by the chain resistance 
R= 1.2 GO. including the resistor inserted between the Greinacher circuit and the cathode and the 
combined capacitance C = 5.5 nF of the Greinacher circuit. The peak current for discharging from 
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Figure 7. Screenshot of the graphical interface of the Greinacher control and monitoring software. 



60 kV is 50 At A. 

The newly developed mechanics have a cable-pull system like a bicycle brake, i.e. the axis rod 
and the feedthrough shaft are connected with two stainless steel wires through two fixed stainless 
steel tubes. Thanks to this new technique the wires can be routed flexibly providing a large freedom 
in choosing the location of the discharging rod independent of the position of the feedthrough. 

Such a discharging system allows not only discharging the circuit but even more importantly a 
direct measurement of the cathode voltage by measuring the peak value of the discharging current. 
The discharging current can be measured by recording the voltage across the 220 kO. reading re- 
sistor inserted between the discharging resistor chain and the ground. The measured peak voltage 
is equal to 0.15 per mil of the cathode voltage, e.g. 1.5 V is measured for 10 kV. A surge arrester 
and a Zener diode are connected in parallel to the reading resistor to protect the measuring device, 
which is usually an oscilloscope. 



3. Characterization of the Greinacher HV multiplier 
3.1 Theoretical model 

The total output voltage expected from a Greinacher chain is ideally equal to the value of the input 
peak-to-peak AC voltage multiplied by the number of stages. Under real conditions however, the 
output voltage can be smaller. Losses are mainly due to effective capacitances of the diodes and 
several stray capacitances which depend on the actual geometrical setup. Altogether the "shunt" 
capacitance is the capacitance between the AC side (the top row in the diagram in Figure ^ and 
the DC side (the bottom row) of the Greinacher circuit. In a transmission line model, which can be 
found in [UPl EOl and also in [1T3I1 the joint effects of charge dissipation in an external load resistance 
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Figure 8. Left: Schematic diagram of the discharging system. Discharging current can be read out by 
measuring the vohage across the 220 k£2 readout resistor Right: Discharging system mechanics. The 
discharge contact and the resistor chain rotates about the polyethylene axis rod steered by the rotary motion 
feedthrough. The axis rod and the feedthrough are connected with two wires through fixed stainless steel 
tubes providing freedom to place the discharging rod independently of the feedthrough position. 



Ri and charge in the shunt capacitances can be calculated to yield the total output voltage, given as 



is the number of stages, v the frequency of the input voltage generator, 7?^ the resistive load (as- 
sumed to be equally distributed on each stage), Cp the shunt capacitance (parallel to each diode) and 
Cs(^> Cp) the capacitance of a single capacitor element. The slope of the Greinacher output equals 
to V^V^pp = tw\\{lyN)/2y. Neglecting the resistive load contribution, the shunt capacitance can 
be estimated directly from it. 

In the actual ArDM setup, where the Greinacher is connected to the field shaping rings and the 
cathode, we can expect surface leakage currents which can affect the response of the Greinacher 
circuit, acting as a Ri load. We first performed a set of measurements in air and in liquid argon 
with the circuit disconnected from the ArDM setup, as described in the next Section. 

3.2 Charging up of the circuit at room and LAr temperatures 

Before assembling the whole TPC we tested the Greinacher circuit alone - not connected to the 
field cage - in air at room temperature and in a polystyrene container open to the atmosphere and 
filled with LAr for characterizing the circuit. For measuring the output HV we used a custom- 
made field mill [ETl . which is a capacitive device allowing measuring the output voltages of the 
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Greinacher circuit without drawing any current. Such a feature is very important since a resistive 
load resuhing in a small current connected to the Greinacher circuit changes its behavior. 

Figure ^ shows thus measured output voltage at the cathode (i.e. the Greinacher stage No. 30) 
while the circuit was charging up. Different curves are for different input AC voltages, i.e. 150 Vpp 
(red) and 250 Vpp (magenta) in air and 150 Vpp (green), 200 Vpp (blue) and 400 Vpp (cyan) in 
LAi\ The horizontal axis shows the time after the AC voltage source was turned on. The input AC 
voltage was ramped up at a fixed rate of 10 V/s throughout the operation of the system. The trend 
of the increase of the voltage at the first part of each curve thus is determined by the ramp-up rate 
and is independent of the input voltage, as can be seen clearly in the zoom for the first 65 s shown 
in the right plot in Figure 0. This trend however is different between in air and in LAr due to the 
threshold effect of the diodes in cold (see also Section ^31 ). For each of the two temperatures it was 
parametrized with a polynomial function and is shown with a dotted curve. 




50 100 150 200 250 300 350 400 10 20 30 40 50 60 



Charging time [s] Charging time [s] 

Figure 9. Cathode voltage recorded as a function of time measured in air at room temperature or in LAr 
while the Greinacher circuit was charging. Different curves are for different input AC voltages, i.e. 150 Vpp 
(red) and 250 Vpp (magenta) in air and 150 Vpp (green), 200 Vpp (blue) and 400 Vpp (cyan) in LAr. The 
right plot is a zoom for the first 65 s. The dotted curves show the trends while the amplitude of the input AC 
voltage was ramping up (10 V/s). The solid curves show the exponential function obtained in Figure [T^- 



To see a general trend of charging of the circuit discarding the effect of the ramp-up rate of 
the AC power supply, the increase of the output voltage after the set value of the input voltage was 
reached is plotted in Figure |Ty. The increase of the voltage is normalized to the difference between 
the voltage measured when the set value was reached and the finally reached voltage. Thus all 
the curves shown in Figure [T^ lie on a single curve and fit very well with an exponential function, 
implying that the charge-up time was approximately independent of the input/output voltages and 
also of the temperature. Fitting all the points with an exponential function V{t)/Vfin = l—e^'i the 
time constant T for the charge up of the circuit was found to be 26 s, i.e. the circuit was charged up 
to 99% of the final value in 120 s after the input AC voltage reached the set value. 
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Figure 10. Increase of the cathode voltage after the amplitude of the input AC voltage reached the set value. 
Different colors correspond to those in Figure 0. All the points agreed with a single exponential function 
shown with the solid line regardless of the input voltage or the temperature. The time constant was found to 
be ~26 s. 

3.3 Linearity of the DC output voltage as a function of the AC input voltage 

Figure [TT] shows the output DC voltage measured at the last stage of amplification as a function of 
the peak-to-peak input AC voltage in air at room temperature (red) and in LAr (blue). The points 
measured in air agreed well with a linear dependence between the output and input voltages over 
the whole range of the tested output voltages up to 7.5 kV. The curve for LAr was shifted towards 
higher input voltage by 33 Vpp showing the threshold effect of the diodes in cold, i.e. the diode 
drew forward cuiTcnt only when the voltage across the three series diodes exceeded the threshold 
value. However a linear dependence was verified also in LAr up to 13.5 kV output although at 
lower voltages the points showed a small deviation, and therefore the points below 100 Vpp were 
excluded from the fit. The slopes obtained from a linear fit were 28.9 VA'^ in air and 28.2 V/V in 
LAr, respectively. In the ideal case the slope equals to the number of the Greinacher stages, i.e. 
30. The shunt capacitances obtained from the slopes (Eq. PTT) ) were 2.7 and 4.4 pF for in air and 
in LAr, respectively. For an identical geometrical configuration, a larger value in LAr is expected 
because of its relative dielectric constant of 1.5. The remaining difference is due to a dissimilar 
capacitance between the Greinacher circuit and the surrounding ground in the two setups (air and 
LAr configuration). 

3.4 Linearity of the voltage distribution across the multiplying stages 

The linearity of the voltage distribution over the field shapers is the first indicator of the uniformity 
of the drift field. Potentials at different Greinacher stages charged at 250 Vpp were measured in air 
at room temperature. A linear behavior of the Greinacher stage voltages was obtained, as shown 
in the left plot in Figure [T^ On the right plot, red points show the measured voltage differences 
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Figure 11. Output DC voltage as a function of the peak-to-peak input AC voltage measured in air (red) 
and immersed in LAr (blue). Same measured in the ArDM setup using the discharging rod (black). See 
Section O. 

between neighboring Greinacher stages. The error shown in the plot corresponds to 0.1% of each 
measured voltage of field shapers. The blue curve is calculated using a transmission line model 
with a shunt capacitance of 2.7 pF as described in the previous section. The expected reduction of 
the voltage across the last Greinacher stage with respect to that of the first stage is ~6%. The mea- 
sured values were in good agreement with expectation although a few points deviated by ~10%. 
The expected reduction can be as large as ~10% in LAr from the calculation based on the shunt 
capacitance of 4.4 pF (Section |?3|) . 

3.5 Characterization of the HV system with the discharging method 

The cathode voltage, measured by discharging the Greinacher circuit using the discharging system, 
as described in Section [Z.3.3[ , can be used effectively to characterize the HV system. A series of 
dedicated measurements were performed with the fully assembled setup in ArDM. 

A typical example of a measured curve is shown in Figure |T3|. The sharp peak at ? = s shows 
the fast discharging of the cathode grid. It is followed by a slower exponential decrease of the 
voltage due to dischai^ging of the Greinacher circuit through the current limiting resistor (200 MQ.) 
in addition to the discharging resistor chain (1 GH), as described in Section |2.3.3| . The discharging 
curve measured across the reading resistor was fitted with an exponential function V{t) = Vb • 
for t > 0.4 s. The intercept Vq is proportional to the cathode voltage before discharging. The shown 
curve was obtained for the circuit charged at 1400 Vpp input AC voltage and Vq = 5.31 ±0. 13 V 
corresponds to the cathode voltage of 35.3 it 0.9 kV. The measured decay constant T = 7.2 it 0.3 s 
is consistent with the naive calculation of T = RC = 6.6 s (see Section [2.3. 3D , taking into account 
possible temperature variations of the capacitors and additional stray-capacitance of the ArDM 
setup. 
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Figure 12. Left: Measured voltage at each Greinacher stage charged at 250 Vpp in air at room temperature; 
Right; Measured voltage differences across each stage and calculated curve (blue) using the transmission 
line model with a shunt capacitance. See text. 
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Figure 13. Discharging curve at the input AC voltage 1400 Vpp fitted with an exponential decrease of the 
voltage. The intercept of the exponential Vb = 5.31 ±0.13 V at f = corresponds to the cathode voltage of 
35.3 ± 0.9 kV. 



To measure the output voltage as a function of peak-to-peak input AC voltage the Greinacher 
circuit was charged up at different input voltages and then discharged recording the discharging 
curve. The measured points are plotted in black in Figure [TT]. Voltage measurements using the 
discharging current clearly had a larger error compared to those using the field mill. The errors 
were dominated by those on the resistances and on the determination of the intercept by the fit and 
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were found to be ~3% of the measured voltages. Within the measurement precision a linearity of 
the output voltage to the input voltage was verified up to the input AC voltage of 1400 Vpp which 
resulted in the output voltage of 35.6 ±0.7 kV (from the linear fit for the input voltages >200 Vpp) 
and the drift electric field of ssO.6 kV/cm. The slope of the linear fit was found to be 26.5 VA'^, 
somewhat smaller than that obtained from the measurements with the Greinacher circuit alone in 
LAr (blue). The observed reduction of the slope by ~ 6% cannot be attributed to a resistive load, 
which is constrained by the observed discharging time (see Section P"!] ) to be Rt > 2.5 TQ.. It can 
be described by two effects: (1) the presence of the DC-voltage shifting circuit which decreases the 
slope by ^3%; (2) a capacitive contribution of the dewai^ inner surface and the field shapers to the 
shunt capacitance. It should be noted that the measured output voltages will be confirmed by the 
drift velocity measurements, as shown in Section 0. 

We also investigated the charging curves at the input AC voltage of 500 Vpp and 900 Vpp. Each 
point was measured in such a way that the input AC power supply was turned off at a certain time, 
while the Greinacher circuit was still charging, and then the circuit was discharged and the cathode 
voltage was measured from the dischai^ging cun^ent. The behavior was found to be similar to that 
shown in Figure ^, with the first part determined by the ramp-up rate of the input AC voltage 
followed by the second part consistent with an exponential function. The time constant of the 
exponential part was also consistent with that obtained from Figure [T^ and showed that the circuit 
was fully charged ~2 minutes after the set value was reached at the input AC power supply. 

3.6 Characterization of tlie HV DC-siiifting circuit 

To check the functionality of the DC-voltage shifting circuit the cathode voltage was measured with 
the discharging method, while the shifting DC HV was applied at different offset values as plotted 
in Figure 0. The measurements agreed with the expected values (dotted curve) verifying a linear 
increase of the cathode voltage depending on the offset voltage within the measurement precisions. 
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Figure 14. Cathode voltage measured using the discharging rod as a function of the DC offset voltage. The 
error bars correspond to ±3% of the measured output voltage. The dotted curve shows the calculated value. 
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4. Operation of the LAr LEM-TPC in the double-phase mode 
4.1 Drift velocity based on scintillation light measurements 

The LAr LEM-TPC was operated in the double-phase mode in pure argon for a continuous period 
of about a month, recording cosmic track events (a large fraction were vertically traversing muons) 
with the PMTs and the LEM charge readout system. The drift electric field was created using 
the built-in Greinacher HV system and was maintained essentially throughout the experiment. The 
Greinacher circuit was re-charged periodically. The recording of cosmic events could be done when 
the charging AC power supply was on or off. The natural decay rate of the HV in the system was 
measured several times during the experiment and was found to be consistent with an exponential 
decay with a time constant of the order of 5-10 days. Since the natural decay rate is very slow, 
changes of the drift field were negligible for a few hours of measurements even when the AC power 
was off. 

With the drift of electrons and their extraction across the liquid-vapor interface, signals of 
ionization in LAr can be detected using the PMTs through the S2 secondary scintillation light. 
It should be noted that the LAr LEM-TPC prototype was not fully optimized for light readout. 
The solid angle coverage by the photocathode of the PMT for the S2 light was only ~0.05%, the 
extraction region being located at the distance of ~70 cm from the PMT, but was sufficient to detect 
clear primary (SI) and secondary (S2) signals, and in particular, their time distribution. Typical 
recorded waveforms for the SI and S2 are shown in Figure [T3[ The two histograms correspond to 
the same cosmic event taken with a drift field of 355 V/cm. The S 1 signal is used to determine 
the event time (Tq = /is). The S2 signal consists of single photoelectron pulses distributed over 
several hundred microseconds, corresponding to the different drift times of the recorded hits. The 
red line in the right figure shows a typical threshold value (corresponding to ~0.25 p.e.) used to 
discriminate PMT pulses from the noise. 

Figure |T^ shows the time distribution of PMT pulses for cosmic events triggered with the 
PMT itself, obtained for four different values of the drift field, i.e. 443 V/cm (magenta), 355 V/cm 
(blue), 266 V/cm (green) and 178 V/cm (red). The vertical axis shows the number of PMT pulses 
above the threshold ~0.25 p.e. regardless of their pulse heights. The pulse counts are normalized 
to the highest peak at ? = /is (Tq), which corresponds to the primary scintillation light (SI). In 
this way the visibility of the S2 signals can be enhanced despite the fact that they mostly consist 
of a bunch of single photoelectron pulses distributed over several hundreds of microseconds. It is 
worth mentioning that, to ensure a pulse linearity for a large signal, relatively large capacitances 
were connected to the four last dynode stages. A polypropylene capacitor having 22 nF was used 
for the dynode 9 and 10, and 122 nF for the dynode 11 and 12. No evidence of a gain reduction 
was seen after a very strong SI pulse (see Figure |T3|) with these capacitors. 

The SI peak is followed by the S2 signal continuing until the maximum drift time which 
depends on the drift field, and can be seen as clear edges in Figure |T^. This maximum drift time 
{tmax) corresponds to the time when the ionization electrons generated in the vicinity of the cathode 
grid, i.e. the bottom of the field cage, arrive at the LAr surface and are extracted into GAr. While 
only those muons which cross both, the top and bottom face of the field cage create S2 signals over 
the full drift time (Tq to fmax)> "grazing" muons which cross one or two of the side walls produce 
an S2 distribution having a smaller width. However, the t^^ax edge is clearly visible when enough 
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Figure 15. A waveform of the PMT signal for the primary (SI, left) and secondary (S2, right) scintillation 
light of the same cosmic event recorded under the drift electric field of 355 V/cm. The red line in the right 
figure shows a typical threshold value (corresponding to ^0.25 p.e.) used to discriminate PMT pulses from 
the noise. The maximum drift time for this event is about 440 /is. 
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Figure 16. Time distribution of PMT signals with a presence of the extraction of electrons across the 
liquid-vapor interface for four different electric drift fields, i.e. 443 V/cm (magenta), 355 V/cm (blue), 266 
V/cm (green) and 178 V/cm (red). The sharp peak at f = /is is the primary scintillation (SI) followed by 
continuous S2 signal. The S2 signal shows a clear edge at the maximum drift time (fmax) depending on the 
drift field. The number of pulses above threshold was counted regardless of their pulse heights to enhance 
the S2 signal. The pulse counts are normalized to the S 1 count. 



events are accumulated, because of the fact that there can be no extraction after ?niax- Since the drift 
velocity Vd{E, T) is a function of electric field (£") and of LAr temperature {T), f^ax also depends 
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on those two parameters 

t^UE,T)-To = —^, (4.1) 

where L is the height of the field cage, i.e. L = 600 mm. The electron drift velocities Vd thus 
calculated from the measured fmax are plotted with black dots in Figure [T7| as a function of drift 
electric field. The drift field was calculated based on the peak-to-peak value of the input AC 
voltage using the linear function obtained from the output voltage measured as a function of the 
input voltage, which is described in Section P3t The horizontal error bars correspond to ±3% 
of the calculated cathode voltage and the vertical ones ±10 /^s in the measurement of ?max- We 
compare our measurements with the ones obtained with ICARUS purity monitors at low electric 
fields [E3I1 taken at 89 K. In order to estimate the drift velocity at the temperature of 86.3 K of 
our setup obtained from the measurements using the PLC system (see Section we employed 
the connection Av^ ~ — 1.7%Ar IHM to scale the fifth order polynomial, obtained from a fit to the 
low field measurements [ E3l1 . Our points agree well with the polynomial function, as shown in 
Figure |T7[ 
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Figure 17. Measured drift velocity as a function of the applied drift field. Black points are measurements 
with proportional light and the red point is a measurement using the charge readout. See text. The dotted 
curve is polynomial function fitted to low field measurements II23II reseated to account for our operating 
temperature. The solid curve shows the WaUcowiak parametrization 1124)1 for E > 0.6 kV /cm computed for 
86.3 K. 



4.2 Drift velocity and uniformity based on ionizing tracks measurements 

The LEM charge readout system operated successfully in double-phase mode. Three examples of 
cosmic events obtained with typical electric field configurations (the drift field of ?»0.4 kV/cm and 
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35 kV/cm in the LEM holes, see Table P are presented in Figure |TE[ Using a simple event display 
images for two projections, i.e. xt- and W -projections, can be obtained for each event: x- and y- 
coordinates are two horizontal coordinates orthogonal to each other and ^-coordinate corresponds to 
the vertical coordinate z- From a visual inspection, the effective gain of the readout is above 15 for 
best field configurations, consistent with results obtained on the smaller detector prototype having 
an active area of 10 x 10 cm^ 0]. A quantitative analysis will be reported elsewhere [ITOl . The 
maximum drift time ?max can be directly determined from those tracks, i.e. the ?max edge is clearly 
visible in a time distribution of the charge signals. The corresponding drift velocity is shown in 
Figure [T7| (with a red dot) and is in a good agreement with the other measured points. 



View Q: Everiit display (run 14450, event BE) 
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Figure 18. Event display for three cosmic events recorded using the LEM charge readout system. The left 
and right columns correspond to the xt- and -projections, respectively. Top: a straight muon track crossing 
the field cage from top to bottom. The time of the start (To ~ 60 /is) and the end (?niax ~ 490 jUs) of the 
extraction are clearly visible. Middle: an example of a "grazing" muon track with delta rays, where the 
extraction starts later at ? ^ 150 jUs. Bottom: interacting event. 

To evaluate the uniformity of the drift electric field along the vertical coordinate, cosmic muon 
tracks recorded using the LEM charge readout system were analyzed. Figure [T9| shows a time 
distribution of the hits for 147 fully reconstructed tracks. The number of readout channels giving 
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a charge signal above threshold (i.e. a hit) within each time bin was counted and put into the 
histogram. The data were taken with the typical field configurations (c.f. Table [I]). The muon 
tracks which crossed both of the top and the bottom faces of the field cage were selected by an 
offline analysis. In addition a continuity of the track was required selecting tracks having no more 
than eight readout channels without a hit along the track. 
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Figure 19. Drift time distribution of hits recorded with the charge readout system for 147 tracks of cosmic 
muons. Continuous tracks crossing both the top and the bottom face of the field cage were selected. The data 
were taken with the typical configurations, i.e. 35 kV/cm in the LEM amplifying holes and a drift electric 
field of 0.4 kV/cm. The red line shows a straight line fit (see text). 
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The number of hits in each time bin (dNbu/dt) is proportional, e.g. in the .x-view, to dx/dt and 
can be related to the electron drift velocity Vd 

dNhit ^dx _dx dz _dx y^^^^-^^-j 2) 
dt dt dz dt dz 
Assuming that the track is straight, the slope dx/dz is constant and hence dNhn/dt is proportional 
to Vd. The distribution in Figure |T^ can therefore be interpreted to be the drift velocity as a function 
of the vertical coordinate z- The drift velocity is constant when the drift field is uniform, i.e. any 
deviation from a fiat distribution indicates a non-uniformity of the field. 

The excess at ? = /^s is an edge effect which can be explained by the strong extraction 
field near the liquid argon surface. The rest of the distribution is straight and can be fitted with 
a line of a given slope to take into account the slight decrease towards a longer drift time, i.e. 
towards the bottom of the field cage. Fitting from 20 to 420 /is and using the average velocity 
(v) = 600/433 = 1.39 mm/jis, one obtained 

Vdit)=vo-at, (4.3) 

with vo = 1.43 ±0.02 mm//is and a = {2.1 ±0.7) • 10""^ mm/jUs^. The drift velocity at the bottom 
of the field cage thus was calculated to be Vd(fmax) = 1.34 ±0.04 mm//is which is 94% of vq. 



From the obtained drift velocity one can derive the drift electric field using our polynomial 
function shown in Figure |r7| (dotted curve). The electric field thus calculated was 392 it 10 V/cm 
at the top of the field cage, respectively 350 ± 18 V/cm at the bottom. The field at the bottom 
was found to be (1 1 it 5)% smaller with respect to that at the top. For the first approximation this 
reduction of the electric field at the bottom of the field cage can be explained by the the non-linearity 
of the voltage distribution to the field shapers due to the shunt capacitance in the Greinacher circuit, 
as described in Section The obtained reduction was found compatible with that of ~10% which 
was expected for the shunt capacitance of 4.4 pF obtained from the measurement in LAr. It should 
be noted that such a non-lineaiity and its influence to the field uniformity can be reduced by further 
optimizations of the design of the Greinacher circuit and of the field cage, respectively. It can also 
be coiTccted for offline. The result reported above demonstrate that this method was effective in 
evaluating the uniformity of the drift field. With increased statistics it will allow a precise correction 
in the conversion of the measured drift time to z-coordinate. 



5. Conclusions 

We have built and operated the first ever LAr LEM-TPC having a large active area with dimensions 
76 X 40 cm^ and a drift length of 60 cm. The chamber was equipped with a 30-stage immersed 
cryogenic Greinacher HV multiplier for generating a HV for the drift field. The characteristics 
of the HV system were extensively studied in air at room temperature and in LAr. An accept- 
able linearity of the output DC voltage to the input AC voltage and of the voltage distribution 
over the Greinacher stages was obtained from the measurements. The LAr LEM-TPC prototype 
operated successfully in the double-phase (liquid-vapor) mode in pure argon during a test period 
of ~1 month, recording three-dimensional images of very high visual quality from cosmic muon 
tracks. The HV system operated stably in the double-phase operation mode of the whole TPC in 
pure argon. As expected the system did not introduce noise problems on the readout electronics of 
the charge readout system, which is one of the advantages of this novel technique promising for a 
giant-scale LAr-TPC. The cathode voltage of w36 kV was reached leading to a drift field of wO.6 
V/cm. 
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